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1 Introduction 

I first got to know of Joe (tfie speaker was MMN) because of his paper with Singh on the 
time-energy uncertainty relation [1]. [See the Appendix for this story.] To use the lan- 
guage of the discussions of this conference, the Eberly-Singh paper was "fundamental" 
though maybe not "useful." But it certainly was fun. That is what we are presenting 
today. Something that is fun. We are going to discuss the Paul trap from a different 
point of view. We will discuss the Paul trap's symmetries and use them to obtain the 
time-dependent solutions, without having to solve the mixed, second-order, partial, dif- 
ferential Schrodinger equation. We will be exemplifying something that Christopher 
Gerry alluded to yesterday: although the Heisenberg and Schrodinger formulations of 
quantum mechanics are equivalent, different things can be more transparent in one of 
them. It was Pauli who taught us this. 



2 The Paul trap 

The Paul trap provides a dynamically stable environment for charged particles [2]- [4] 
and has been widely used in fields from from quantum optics to particle physics. 

Paul gives a delightful mechanical analogy [4] . Think of a mechanical ball put at the 
center of a saddle surface. With no motion of the surface, it will fall off of the saddle. 
However, if the saddle surface is rotated with an appropriate frequency about the axis 
normal to the surface at the inflection point, the particle will be stably confined. The 
particle is oscillatory about the origin in both the x and y directions. But it's oscillation 
in the z direction is restricted to be bounded from below by some zq > 0. 

The potential energy can be parametrized as [2] 

V{x, y, z, t) = 14 (z, t) + Vy{y, t) + V, (z, t), (1) 

where 

V^x.t) = +-l^V{t)x^ = g{t)x^, (2) 



2r? 

e 







yyiy^t) = +^V{t)y'^g{t)y\ (3) 

V,iz,t) = -4 V(i) ^' = .93(0 (4) 

These potentials can be used to set up the classical motion problem 

x,i=F,{x,t) = -^YlpH = ^2g{t)x,i, (5) 

ax 

y.i^Fy{y,t) = -'^^y^ = ^2g{t)y,u (6) 

z,i = F,{z,t) = -^^-^^ = -2gs{t)z,i. (7) 

The solutions to these equations are oscillatory Mathicu functions for the bound case 
[2]. The oscillatory motion is similar in the {x,y) directions but is different in the z 
direction, since g{t) ^ gait). 



3 The quantum-mechanical Paul trap 

In the quantum mechanical treatment of the Paul trap, one must solve the Schrodinger 
equation, which in Cartesian coordinates has the form (h = m = 1) 

S^{x, y, z, t) = {a,, + dyy + d,, + 2idt - 2g{t){x^ + y^) _ 2g^(t)z^] ^[x, y, z, t) = 0. 

(8) 

The time-dependent functions g and .93 in Eq. (8) are 

9{t) = + A vw, 53W = -4 nt)- (9) 

where 

V{t)^Vdc-VacCOSUj{t-t^) (10) 

is the "dc" plus "ac time-dependent" electric potential that is applied between the ring 
and the end caps. 

Exact solutions for the 1-dimensional, quantum case were first investigated in de- 
tail by Combescure [5]. In general, work has concentrated on the z coordinate, but not 
entirely [6] . Elsewhere [7] , stimulated by the work of Ref. [8] , we discussed the classi- 
cal/quantum theory of the Paul trap in the z-dircction. 

The form of the Schrodinger equation (8) suggests that, in addition to the above 
Cartesian coordinate system, there is another natural coordinate system for this prob- 
lem: the cylindrical coordinate system. Introducing such a change of variables, 

a; = rcos6', y = rsiiiO, z = z, (11) 

the Schrodinger equation becomes 

Scyi<^{r, e, z, t) = \drr + '^dr + ^dge + 8,, + 2idt - 2g{ty - 2g3{t)z^^ $(r, 9, z, t) 
= 0. (12) 

In addition, the forms of Eqs. (8) and (12) suggest that we can factorize the solutions 
and equations into the forms 

^{x, y, z, t) = X{x, t)Y{y, t)Z{z, t), (13) 

S,X{x,t) = {d.,^ + 2tdt-2g{t)x^}X{x,t)^0, (14) 

SyY{y,t) = {dyy + 2idt " 2g{t)y^} Y{y,t) = 0, (15) 

S,Z{z, t) = {d,, + 2idt - 253(t)^'} Z{z, t) = 0, (16) 

and 

$(x,y,z,t) = 17(r,0,t)Z(z,t), (17) 

Sren{r,e,t) = \drr + -dr + \dee-2g{ty\n{r,e,t) =0. (18) 
y r r'^ J 

Being physicists, we tend to blindly go ahead and accept this separation, ignoring the 
mathematical subtleties in separating coordinates in time-dependent partial differential 
equations. Fortunately, Rod insists on keeping me honest. But this procedure does turn 
out to be justified [9]. So, now we can just go blindly on. 



4 Lie symmetries and separable coordinates 

Lie symmetries for the Sclirodinger equation (8) can be obtained by solving tlie operator 
equation [10, 11] 

[S,L]^X{x,y,z,t)S. (19) 

The operator S is one of the the Schrodinger operators we have discussed, L is a gener- 
ator of Lie symmetries, and A is a function of the coordinates x, y, z, and t. An operator 
L has the general form 

L = Codt + CA + C2dy + C^d, + C, (20) 

where the coefficient in each term is a function of the coordinates and time. 
First we define what are going to be useful separable coordinates: 



y 



t, (21) 



and 



2 



The t-dependent functions (f>(t) and 03 (i) are given by 

0==2ee 03 = 26f3, (23) 

where {^(t), £,it)} and {^3(0: ^3(^)} the complex solutions of the second-order, 
linear, differential equations in time 

7 + 2(7(^)7 = 0, 73-^253(^)73=0, (24) 

respectively, and satisfy the Wronskians 

W{C0=W{^s,^3)^-i. (25) 

Equations (24) are the same as the classical equations of motion (5)- (7) for the Paul 
trap. The solutions are the same Mathieu functions. For certain values of the parameters 
in the potential (1), the solutions are bound, meaning the charged particle is indeed 
"trapped" [2]. This shows a connection between classical and quantum dynamics. 

5 Cartesian symmetries 

For this exercise we can concentrate only on the z coordinate, since formally the results 
are the same for the x and y solutions, with the exception that the ^3's and gs's, etc., 
lose the subscripts 3. [Elsewhere, we will discuss the symmetries of the Paul trap in 
much greater detail [9].] 

One can find, or simply verify, that the Schrodinger equation has the symmetry 
operators 

Jz- = - = +^ (I) ' [az + z (1 - 103)] , (26) 
J.+ = -hd. + *e3^ = (I) ^ [az - z (1 + §03)] . (27) 
These operators satisfy the nonzero commutation relation 

[J,_,J,+]=/, (28) 



and so form a complex Heisenberg Weyl algebra w'^. This means that the operators 
generate a set of "number states" given by 



(29) 
(30) 



There is also an su(l, 1) algebra, which we will not go into here. This is a generalization 
of the "squeeze algebra." 

A word of caution. The states Zn, {z, t) are not to be construed as energy eigenstates. 
Zn^ is a solution to the time-dependent Schrodinger equation (16). It is generally not 
an eigenfunction of the Hamiltonian. That is, 



idtZn^ = HZn, = i 



Zn^ 



^ [Const] Z„^. 



(31) 



So, to keep John Klauder from getting mad at me, we use the terminology "extremal 
state" for Z^ and "higher-order" states for Z„, . We restrict the terms "ground state" 
and "excited state" for problems which allow eigenstates of the Hamiltonian. 

Eq. (29) is a simple first order differential equation for Zq: 



J,-Zo{z,t) = 0. 



Solving it yields 



Zo(z,i) = /(t)cxp. 



1 



(32) 



(33) 



where the proportionality constant, / must be a function of t. 

A first idea might be to take this proportionality constant as the "normalization 
constant," [n(l)3{t)]^^/'^ . This would conserve the probability, as one would want for the 
time- development of a unitary Hamiltonian. But this turns out to be incorrect. Such a 
solution does not satisfy the Schrodinger equation (16). Indeed, putting Eq. (33) into 
Eq. (16) yields a first order differential equation in t for f(t). The resulting normalized 
extremal-state solution is 



Zo{z,t) = (TTi^ia)" 



-1/4 



a/4 



exp 



(34) 



So why is the phase factor (63/^3) there? It is there because, in this time- 
dependent Schrodinger equation, it is the necessary generalization of the simple- 
harmonic oscillator ground-state energy exponential, eKjp[—ihio/2]. (This follows since 
^3(0 (2tj)~^/^ exp[icji].) This phase factor is necessary [12, 13, 14] for Eq. (34) to be 
a solution of Eq. (16). 

Now repeated application of Eq. (30) gives 



Z„^iz,t) = [nj]-'/^ [J,+r^ Zoiz.t). 



(35) 



But the right hand side of Eq. (35) is proportional to [9] the Rodrigues formula for the 
Hermitc polynomials [15]. Using this yields the result 



Z„^(z,t) = (2"=nJ)-^/'(^03)"'^' 
i/„^(z) exp 



M"^ + 5) 



Z 



(36) 



The forms of (x, t) and Y„ {y, t) follow immediately by just changing notation, 



Xn^ (X, t) 



so that 



(37) 
(38) 

(39) 



(a;, y, z, t) = (a;, t)r„,^ (y, t)Z„^ (z, i). 
6 Polar symmetries 

In what is intuitively interesting, finding the symmetries for the polar coordinates 
amounts to considering complex linear combinations of the Cartesian operators. Take 
the following two pairs of operators which form the basis of two Heisenberg-Weyl alge- 
bras [11]: 

]^{J^-+iJy-) = \j\ i{dx + idy) - iS,{x + iy) 



1/2 



dp + -dg+ p[l- 



1/2 



P 



S,{dx - idy) - ii{x - iy) 



If we add to the above four operators /, the operator 

K. = Jx+Jx- + Jy+Jy- + 1 = a+a- + c+c_ + 1, 
and the angular momentum operator 

Cz = i{ydx - xdy) = ~ide, 



(40) 

(41) 

(42) 
(43) 



we have that the set {I,a±,c±,IC,Cz} forms a closed algebra. Further, if wc make the 
transformations 



/ = i (/C - £,) = a+a_ 



1 

2 ' 

d = i (/c + /:^) = c+c_ + i 



(44) 
(45) 



we see that the the symmetry algebra has two oscillator subalgebras, {/, a±,/} and 
{d, c±,/}, which have only the identity operator, /, in common. Therefore, the algebra 
S'x y = osa(l) + osc(l), with the two Casimir operators Ci = a+a- — / = — ^ and 
C2 = c_|_c_ - d = -i. 

We restrict ourselves to the representations of os(l) that are bounded below; namely, 
1-if2 + T-1/2- Let n„^,„ be a member of the set of common eigenfunctions of / and d 
spanning the representation space. Then, for {n,m) G 7iq, we have 





d^n.m ^2) 


(46) 


n— l,m ; 




(47) 


71-^1, 7n 7 




(48) 


(/ - d)nn 


ra (j^ ^)^7i,mi 


(49) 




= (n + 7TI + l)r2„,m. 


(50) 



In order that the spectrum of d and / be bounded below, we also have 

a^Qo = 0, c_51o,o ~ 0. 



(51) 



Eqs. (43) and (49) immediately tell us that the (normalized) 9 dependence of the 
^^n,m(f, d, t) can be given by 



f^n,m,(^. ^':*) = "^n^m (?', i) B„,to (6*) =7?.„,m(r, i) 



exp[i(TO — n) 



(52) 



Therefore, ^o_Q{r,9,t) is a function only of r or p and t. In this {n,m) = (0,0) case, 
both the equations in (51) have the same first-order differential form. So, we can solve 
for f2o,o similarly to as was done for the solution. The normalized solution to the 
Schrodinger equation (18) is foimd to be 



f^o,o('^, 0,t) 



roii/2 



'2ti 



exp 



2 



= Xo(x,t)yo(y,t). (53) 



Repeated application of Eqs. (48) gives us the general result: 
f^n,m(p,e,i) = (n!TO!)-^/'a'|c™17o,o(p,e,i). 



(54) 



This turns out also to be a Rodrigues-type formula, although more complicated than 
before [9]. This time it is for the generalized Laguerre polynomials [15]. The end result 
turns out to be 



exp[i(m - n)6l] (-)'^fc! / 2 V'"^ ( 1 



pl"-'"l4l"-™l'(p2)cxp 
i(n + m — |m — n|). 



1 - i 



^(n+m+l) 



2 



Now let us change to the more standard cylindrical quantum numbers: 
€ = 1^2 1, £z ~ m — Ur ~ m + n > 0. 

We find that 

e,,(0) 

V ^TT 



cxp[i £z 0] 



V27r 

(_l)("^-«)/2 



(55) 
(56) 

(57) 

(58) 
(59) 



1/2 , i(„+m+i) 



P ^1 



(60) 



Except for our minus-sign phase convention, Eq. (60) resembles the standard result for 
the ordinary two-dimensional harmonic oscillator [16]. 

This standard solution has the known property that and 1^ cannot differ by an 
odd integer to allow a normalizable solution, [^(n^ — €) must be an integer.] If one 



solves the Schrodinger equation directly, this falls out, but the physical reason is not 
transparent. However, from the symmetry point of view, the reason is clear. The n and 
m quantum numbers, which reflect the fundamental symmetries, can be all non-negative 
integers. The Uj. and quantum numbers reflect a rotation of the axes by 45°. (See 
figure 1.) So, reaching the allowed positions of quantum numbers along these diagonals, 
scaled by l/-\/2, means all integer values of {nr,£z) are not allowed. 




Fig. 1. A plot of the allowed quantum numbers for the polar coordinates of the Paul 
trap. Shown are the (n, m) quantum numbers of Eqs. (54) and (55) as well as the 
{nr,iz) quantum numbers of Eqs. (57) and (60). 



This is similar in spirit to Pauli's being able to show, by symmetry methods, that 
there is an extra conserved quantity, the Runge-Lenz vector, om the hydrogen atom. 
This allowed an understanding of the extra degenerate quantum number. 

Finally, we can write out the complete solution in spherical coordinates as 

*«.A,«Jx,2/,z,i) = i?„„,f(r,i)efJ0)Z„^(z,t). (61) 
Appendix: Getting to know (about) Joe 

Peter Carruthers, who was my advisor at Cornell, had a great ability to find talent in 
physicists. He also was very perceptive about the qualities of people, being a well-known 
lover of life. For reasons that will become clear below, let me give you two examples of 
Pete's first ability. 



i) Wc graduate students were amused by this new. young, hot-shot, assistant pro- 
fessor who was going to pubhsh a famous paper. Wc were more amused when he made 
associate professor. Later, I heard of the move to make him full professor. (From lack 
of being there some of these details must be off.) It was boiling down to Hans Bethe's 
call. Pete figuratively pounded Hans' desk top telling Bethe he had to promote this 
guy. In the end the youngster got his full professorship, published his paper (on the 
renormalization group), and the rest is the history of Ken Wilson. 

ii) When Pete became division leader at Los Alamos. Pete decided to hire as a staff 
member a young guy who had had an unproductive postdoc career - solely on the basis 
that Pete sensed creativity. Pete had to hide him in the Division Office, not a group, 
because funding could not be justified. He just sat around the corner, playing with his 
HP calculator, getting funny numbers again and again, leading to the period-doubling, 
chaos revolution of Mitchell Feigenbaum. 

I mention this all in the context of when I first heard of Joe Eberly. (It was not when 
I was at Cornell, because then my main contact with Rochester was with the wild man 
of the time, Robert Marshak.) It was when Pete came out to Los Alamos. I asked him 
if he had seen the new paper by Eberly and Singh [1] on the time-energy uncertainty 
relation that had referred to our work. Pete said he had, and had liked it. Then he 
remarked that Eberly "is very smart." Being properly impressed, I figuratively wrote 
down Joe's name in my black book of smart guys. Then Pete said he's also a nice a guy. 
So, Joe also went into the "nice guy" book. Most of us, if we are honest, must admit to 
really wanting to be in the first book. But Joe is in both books, and that is why we are 
all here today to honor him. 
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